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We present the asymptotic distribution theory for a class of increment- 
based estimators of the fractal dimension of a random field of the form 
g{X{t)}, where (; : R ^ R is an unknown smooth function and X{t) 
is a real-valued stationary Gaussian field on R'', d = 1 or 2, whose 
P^ ■ covariance function obeys a power law at the origin. The relevant the- 

C/j ' oretical framework here is "fixed domain" (or "infiU") asymptotics. 

(-H , Surprisingly, the limit theory in this non-Gaussian case is somewhat 

richer than in the Gaussian case (the latter is recovered when g is 
affine), in part because estimators of the type considered may have 
an asymptotic variance which is random in the limit. Broadly, when 
g is smooth and nonaffine, three types of limit distributions can arise, 
types (i), (ii) and (iii), say. Each type can be represented as a random 
^ ' integral. More specifically, type (i) can be represented as the integral 

ly^ I of a certain random function with respect to Lebesgue measure; type 

^S) . (ii) can be represented as the integral of a second random function 

lO ' with respect to an independent Gaussian random measure; and type 

^^D , (iii) can be represented as a Wiener-Ito integral of order 2. Which 

>— P ■ type occurs depends on a combination of the following factors: the 

roughness of X{t), whether d = 1 or d = 2 and the order of the incre- 
ment which is used. Another notable feature of our results is that, 
f^ ' even though the estimators we consider are based on a variogram, 

"t^ I no moment conditions are required on the observed field g{X{t)} for 

a the limit theory to hold. The results of a numerical study are also 

presented. 

> 

/\^ • 1. Introduction. 

1.1. Background. The problem of quantifying the roughness of a (con- 
tinuous but rough) curve or surface, whose height is observed at discrete 
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locations on a rectangular grid, arises in many areas of science and technol- 
ogy. A widely used approach is to model the curve or surface as a random 
field whose covariance function follows some form of power law behavior at 
the origin, and then to estimate a scale-invariant measure such as the frac- 
tal dimension to quantify roughness. A good entry point to recent statistical 
literature on this topic is the discussion paper by Davies and Hall (1999). 

Frequently in this approach it has been assumed that the curve or surface 
is Gaussian. However, in a number of applications the Gaussian assumption 
may be open to doubt, and therefore the problem of estimating the fractal 
dimension in non-Gaussian settings is of interest. In this paper we study 
estimators of fractal dimension for a class of stationary non-Gaussian ran- 
dom fields, and we provide a detailed account of the asymptotic distribution 
theory for these estimators, as well as studying their numerical properties 
using simulation. As a preliminary, we provide a brief review of recent work 
in which the Gaussian assumption is made. 

1.2. The stationary Gaussian model. For simplicity we focus mainly on 
the one-dimensional case (d = 1) in this Introduction. Let {X{t):t G R} 
denote a stationary Gaussian process with a covariance function 7 which 
obeys the following power law at the origin: 

(1.1) 7(t)=7(0)-c|t|° + O(|t|"+^) as|tHO, 

where a £ (0,2], known as the fractal index, governs the roughness of the 
sample functions and is typically the parameter of greatest interest in rough- 
ness studies; the positive quantity c is a (local) scale parameter known as 
the topothesy; and /3 > governs the size of the remainder term in (1.1). 
There is a simple relationship, under model (1.1), between a and the fractal 
dimension D of the graph of the sample function, given hy D = 2 — a/2; see, 
for example, Adler (1981). [This result generalizes to D = d + 1 — a/2 when 
X{t) is a stationary Gaussian random field on R"^ with covariance function 
of the form (2.3).] Thus, the larger the value of a, the smoother the sample 
function. 

Suppose we observe a sample 

(1.2) Sn = {X{i/n) -A = 0,1,..., n-1} 

of observations of X{t) at equally spaced locations in the region [0,1]. In 
this formulation of the problem the asymptotic regime as n — > 00 is known 
as "fixed domain asymptotics" and is often appropriate when interest is 
focused on roughness at fine scales. See Stein (1999) for further discussion 
of this type of asymptotic regime. 

A variety of estimators of a, based on data 5„,, have been studied in recent 
years under model (1.1), with X assumed to be stationary and Gaussian. 
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For example, Hall and Wood (1993) considered box-counting estimators; 
Jakeman and Jordan (1990) and Constantine and Hall (1994) discussed 
estimators based on the variogram; Feuerverger, Hall and Wood (1994) 
considered estimators based on counting upcrossings; and Chan, Hall and 
Poskitt (1995) considered estimators based on the periodogram. Kent and 
Wood (1997) considered two modifications to variogram-based estimators: 
the use of higher-order increments and the use of generalized least squares. 
The use of higher-order increments in a closely related context was investi- 
gated independently by Istas and Lang (1997). 

In this paper the focus is on variogram-based estimators of fractal dimen- 
sion, and we now discuss these in more detail. The theoretical variogram 
is given by u{h) = E{X{t + h) - X{t)Y and under (1.1) u{h) ~ 2c|/i|" as 
\h\ — > 0. Given data of the form (1.2), we can estimate u{u/n) by 

(1.3) l>„ = ^[X{(^ + u)/n}-X(^/n)]^ 

i 

and we can estimate a using log-log regression based on the approximation 
relationship 

(1.4) logUu^ const. + a\ogu. 

The simplest estimator of this type is the ordinary least squares estimator 
given by 

,, ,. . ^ E:r=i log ^u (log u-m-^ Er=i log fc) 

^■^ " E™=i(logn-m-iEr=ilogA:)2 ' 

where m stays fixed as n — > oo . 

It turns out that when < a < 3/2 the estimator (1.5) has variance of 
order n~^, but when 3/2 < a < 2 the variance of (1.5) is of order n^"~^. See 
Constantine and Hall (1994); related results were also obtained by Jakeman 
and Jordan (1990). More recently, it was noted by Kent and Wood (1997) 
and Istas and Lang (1997) that if one bases the variogram in (1.3) on second- 
order differences, that is, 

(1.6) Uu = Y.[X{{i + u)/n} +X{{i- u)/n} - 2X{i/n)]', 

i 

then the variance of the resulting estimator (1.5) is of order n~^ for all 
aG(0,2). Thus, there is motivation for considering variograms based on 
increments (i.e., differences) of higher order. In our terminology (1.3) is 
based on an increment of order 0, and (1.6) is based on an increment of 
order 1, and in either case (1.3) may be written in the form 

2 
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where the a^ notation for increments is described in detail in Kent and Wood 
(1997) and Chan and Wood (2000) and is summarized in the Appendix. 

More recently, the corresponding estimation problem for the two-dimensional 
case was considered in the discussion paper by Davies and Hall (1999). The 
analysis of data obtained by two-dimensional sampling is trickier than the 
one-dimensional case, in part because of the possibility that anisotropy is 
present. A helpful discussion of this issue is given by Davies and Hall (1999). 
A second approach to the analysis of two-dimensional surface data is con- 
sidered by Chan and Wood (2000). 

1.3. Non-Gaussian data. In this paper we assume that, instead of ob- 
serving the stationary Gaussian process X{t), we observe a stationary non- 
Gaussian process g{X{t)}, where g is a smooth but unknown nonafiine func- 
tion and X is a stationary Gaussian process satisfying (1.1) as before. Two- 
dimensional sampling of random fields (corresponding to d = 2) is also con- 
sidered. We address the following question: what are the asymptotic prop- 
erties of the variogram-based estimators when we observe g{X{t)} rather 
than X{t)l Prom a theoretical point of view, the quantity being estimated, 
namely the fractal dimension of the sample function, will be the same as 
that of the underlying Gaussian random field provided g is nondegenerate 
and reasonably smooth; see Hall and Roy (1994) for the relevant results. 
However, it turns out that the asymptotic distribution theory for nonaffine 
g is somewhat richer than in the Gaussian case (though the Gaussian case 
is recovered, of course, when g is an affine transformation). 

Chan and Wood (2000) obtained (correct) preliminary results concerning 
rates of convergence in the case when g{x) = x^. However, no concrete results 
concerning the asymptotic distribution theory were given, and in fact the 
conjectures in Remarks 5.5 and 5.6 of that paper do not adequately describe 
the results given below. Also, it turns out that the case g{x) = x^ is not fully 
representative of general smooth nondegenerate g, because some components 
of the limit distribution disappear or are constant in the former case. 

Our main theoretical results are stated in Section 2. One general point 
which emerges is that the estimation of fractal dimension in this non-Gaussian 
setting is, in a sense, more difficult than in the Gaussian case; see, in partic- 
ular, point 2 in Section 2.3. However, on the numerical side a fairly extensive 
simulation study reported in Section 3 suggests that, in practice, the dete- 
rioration in the non-Gaussian setting is fairly mild, provided that g is not 
drastically nonaffine in the relevant domain. The main theorems are proved 
in Section 5. These proofs make use of several lemmas which are proved in 
Section 4. 

2. Main results. The principal results in the paper are presented in The- 
orems A and B. Theorem A covers those cases in which A + Ap—2a>d and 
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Theorem B covers those cases in which 4 + ip — 2a < d, where p > is the 
order of the increment used (see the Appendix), d= 1,2 is the dimension of 
the parameter set of the underlying random field and a G (0, 2) is the fractal 
index. 

Throughout this paper we make frequent use of the notation for incre- 
ments and multi-indices given in the Appendix. 

2.1. Preliminaries. Let {X(t) : G R } denote a real-valued stationary 
Gaussian process {d = 1) or field {d = 2) with covariance function given 
by 

^{t)=cov{X{s),X{s + t)}. 

Let 5 : R ^ R be an unknown function. It is assumed that we observe 
g{X{t)} rather than X{t). 
Define the index set 

(2.1) 2:„ = {0<j<no}. 

When d = 1, j and no are integers and we take n = hq; and when d = 2, 
j and riQ are multi- indices in Z^ (see the Appendix) and no = no(n) is a 
sequence such that n is the product of the elements of no, that is, n = 
no [1] no [2]. Thus, n is the number of elements in In- The dataset we actually 
sample is given by 

■Sn = {gi = g{X{i/no)} : -mJ < i < no + mJ}, 

where division of multi-indices j, A; E Z*^ is defined by j'/A: = (j[l]/A;[l], . . . ,j[d]/k[d]), 
assuming k[l] 7^ for each I. In the definition of Sn, J is a multi-index which 
depends on the increment that we use, and m, an integer, is the number of 
dilations of the increment that we consider [see the Appendix and also Kent 
and Wood (1997)]. It is assumed throughout that m stays fixed as n — > 00; 
justification for keeping m fixed is given by Constantine and Hall (1994). 
Note that the sampling regime indicated above corresponds to "fixed do- 
main" asymptotics as n — > 00. 

Consider the following conditions on 7 and g. 

{Al)q For some a G (0, 2) and /? > 0, and for each nonnegative multi-index 
r with \r\ = a, 



(2.2) 



i^'^Ht) = -_{||trM(t/iit||)} + omr"-^-') 



as \\t\ 



where ||t|| = (t^t)^'^ is the usual Euclidean norm on R'^, and, for a 
nonnegative multi-index r = (r[l], . . . ,T[d]), "y^^'{t) = d'^''y[t)/dt^ = 
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d\'^h{t)/dtf^^---df/\ where \r\ =Ef=i^B- In (2-2) M(-) is as- 
sumed to be a positive constant when d = 1; when d = 2 M(-) is 
a function on the unit circle in R^ , all of whose partial derivatives 
derivatives of order q+1 are assumed continuous. 

(^2) The seventh derivative of g, g^"^', is continuous on R. 

(^3) The set {x:g^^'{x) =0} has Lebesgue measure 0, where g'^^' is the 
derivative of g. 

(AA) {d = 2 only.) As n — > oo no[l]/no[2] stays bounded away from 
and oo. 

It is easily shown that, if (^l)g holds, then (.41)}. holds for each 1 < 
r < q. Moreover, {Al)q implies that 

(2.3) 7(i)=7(0)-||t|rM(i/||t||) + O(||t|r^), 

where (2.3) reduces to (1.1) when d=l. 

It is possible to weaken assumption (^2) to some extent, but only at a 
considerable cost in technical detail in the proofs. Assumption (.43) is a mild 
nondegeneracy condition which seems essential if our theorems are to hold. 
Condition (.44) is needed when d = 2 to ensure that the sampling set does 
not become too "thin." 

Putting gi = g{X{i/nQ)}, we define 

(2.4) Z„ = n-i^|^aj5.+,| , /i„ = n-/'^i?|^ajx(^^) | 

and 

(2.5) Gr= f [g^'Hxmf'dt, r = l,2, 

where d = 1 or 2. The notation for increments a = {oj} used in (2.4) is 
explained in the Appendix. Let L„, u = 1, . . . ,m, be real numbers which 
satisfy 

m m 

u=l u=l 

Various choices for the L^ are discussed by Kent and Wood (1997). All 
are based on the log-linearity given by the power law relationship; see (2.3) 
and (1.4). The simplest case is ordinary least squares [see (1.5)] for which 

(logn-m-^E"=ilog^) „_. 

E«=i(logu-m ^J2y=ilogvy 

Define 

m m 

(2.6) a = ^Lu log Zu and an = ^Lu log //„ 

u=l u=l 
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and 

(2.7) F=(Fi + F2)/Gi, 



where 



(2i 



and 



^1 = S X! -^«^0,l('^0,l« + 7"0,2« + T0,3«) 



.u=l 



te[o,i]'' 



/ , ^uA^o,m(''"0,1« + ■'"0,2m + T0,4m) 



>«=1 



(2.9) 

X 

In the above /xo,n = hm^^oo /"« G (0, oo) and roj« = hm„^oo Tju (j = 1, . . . , 4; 
u = 1, . . . ,7Ti), where fiu is defined in (2.4) and ti„, r2u, rs^ and t^u are de- 
fined in (5.13), (5.18) and (5.21). Note that /xo,n and tqju may be determined 
exphcitly using (1.1) when d = 1 or (2.3) when d = 2, but we omit the details 
because these formulas are not required in what follows. 

The quantity a which appears in Theorem A is defined as follows. When 
d = l. 



(2.10) cJ = VFi^, 

where tu = C^^u~"Lu, n = 1, . . . , m, t = (ti, . . . , tm)'^ , C > is the constant 
in formula (2.7) in Kent and Wood (1997) and <&o is the covariance matrix 
defined via (3.3) and (3.6) in Kent and Wood (1997) (the precise definitions 
of C and $o need not concern us here). When d = 2, a is still of the form 
(2.10), but with the quantities C and $o now given by (3.13) and (3.16) in 
Chan and Wood (2000). 

The integrator B{-) which appears in Theorem A is a Gaussian random 
measure such that, for any measurable subsets Ai,. . . ,Ak of [0, 1]"', 

(2.11) {B{Ai), . . .,B{Ak)f ~ Nk{^, ^), 

where ^ = {'4'ji): i^ji = ^{^j HA;) and A denotes Lebesgue measure on R . 
In Theorem B Zs is the zero-mean random Gaussian measure defined on 
[— 7r,7r] with the following properties: if Di,D2 ^ [— 7r,7r] , then 

cov{Zs{Di), Zs{D2)} = S{Di n D2), 

where S is the spectral measure of a covariance function on Zi'^ of the 
form p{k) = \\k\\~°'A{k/\\k\\), where A; G Z and A{-) > is a continuous 
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function on the unit sphere in R . In fact, p{k) is defined as the Umit as 
n ^ CO of the covariance cov(yi, Fj+fc), where Yi = Yl^=i ^'^^iu/ iJ2^=i u'^y^'^ 
and Yiu is defined in (5.1). When d=l, we may take yl to be a constant, 
and when d = 2, A is determined by M in (2.3). 

2.2. The theorems. We are now ready to state our main results. Theo- 
rems A and B. Discussion of these results follows in Section 2.3. 

Theorem A. Suppose that 7 and g satisfy conditions (^1)4 , (-^2) 
and (.A3), and, if d = 2, suppose also that the sampling regime satisfies (A,4). 
Let a and an be defined as in (2.6), using the index set (2.1) in (2.4). 
Suppose that a = {oj} is an increment of order p>0 such that A + Ap — 2a> 
d, where d = \ or d = 2. Then 

(2.12) a-an = n~^l'^F + n~^'^a^^Z + Op{n-''/'^ + n'^/^), 

where a > is a constant and Z r^ N{0, 1) is independent of Gi, G2 and F, 
where the latter are defined by (2.5) and (2.7)-(2.9). 

When <2a < d, the n~°^''^ term in (2.12) is dominant and, in this case. 



i"/'^(a-a„)^F. 



When d <2a < min(4,4 + 4p — d), the n ^'^ term in (2.12) is dominant 
and, in this case, 

ny\a-an)^aGi' f [gW {X {t)}f dB (t) ^ a^ Z , 
J[0,1]'' Gi 

where B{t) is the Gaussian random measure described in (2.11). When 2a = 
d, both terms contribute, and we have 

n^''^{a-an)^F + aG7^ I [g^^'^ {X{t)]f dB{t). 

"'[0,1]'' 

In the above the Gaussian measure {B(t)} is independent of the underlying 
Gaussian field {X{t) : t G [0, 1]"^} . 

When 4 + 4p — 2a <d, the limit distribution of a may be expressed as a 
Wiener-Ito integral of order 2. See Dobrushin and Major (1979) and Ma- 
jor (1981) for further details on Wiener-Ito integrals. 

Theorem B. Suppose that 7 and g satisfy conditions (^.1)4, (A2) 
and (A3), and (when d = 2) the sampling regime satisfies (A4). Let a and an 
be defined as in (2.6) where, as before, a= {oj} is an increment of order 
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p = 0, and suppose that either d=l and 3/2 < a <2 or d = 2 and 1 < a <2 
(corresponding to the condition A + Ap — 2a <d). Then 



Jxi,x2e[-TT,n]'i Jte[o,i]'^ 

x[g'^^^{X{t)}fdtdZs{xi)dZs{x2) 

= / I{xi^X2)J^g{xi+X2)dZs{xi)dZsix2), 

where the indicator function I{xi ^ X2) excludes the diagonal and 

Jte[o,i]'* 

is the Fourier transform of the random probability measure on [0, 1] whose 
density with respect to Lebesgue measure is given by [g^^'{X{t)}]'^ /Gi. In the 
above, the Gaussian measure Zg is independent of the underlying Gaussian 
field {X{t):t£[0,lY}. 

2.3. Discussion. A number of comments concerning Theorems A and B 
now follow. 

1. When g is affine in Theorem A, F = 0, \fG2lG\ = 1 and therefore the 
limit distribution of n^''^{a — a„) is A^(0,cr^). This agrees with the cen- 
tral limit theorems given in Kent and Wood (1995, 1997) and Chan and 
Wood (2000). 

2. Note that, in (2.12) \fG2lG\ > 1 by the Cauchy-Schwarz inequality, with 
equality if g is affine. Moreover, when a < d/2 the rate of convergence of 
a to an is of slower order than n~^'^. Therefore, we may conclude that, 
from the point of view of the estimator a in (2.6), the non-Gaussian 
case {g not affine) is less favorable than the Gaussian case (^f afhne) in 
the framework considered in this paper. This finding is confirmed by the 
numerical MSE results in Section 3, though these results also suggest 
that the deterioration is not too severe provided that g is not drastically 
nonaffine over the relevant domain. 

3. At the borderline between Theorems A and B [i.e., when the increment 
used has order and a = (4 — d)/2], the limit distribution is of the type 
given in Theorem A, but the convergence rate of n~^'^ is modified by 
a logarithmic factor in n; compare the Const antine-Hall (1994) result 
when a = 3/2. We omit the proof. 

4. The bias of a in the context of Theorem A depends not only on the right- 
hand side of (2.12), but also on a„ — a. It follows from (2.1) and (2.3) 
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and the definitions in (2.6) that a„ — q = 0{n~"' ). It can be seen that, 
in Theorem A, this term makes a negligible contribution to the bias if 
and only if /3 > a, and it makes a negligible contribution to the mean 
squared error if and only if /? > min(a, d/2). In Theorem B a„ — a makes 
a negligible contribution to the MSE if /3 > 2 — a. See Constantine and 
HaU (1994), Kent and Wood (1997) and Chan and Wood (2000) for fur- 
ther discussion of this bias term. 

5. It is interesting to note that Theorems A and B do not require any mo- 
ment conditions on g{X{t)}. This is a consequence of the Gaussianity of 
X{t) and the smoothness of g. See also the proof of Step 1 in Section 5. 

6. It is straightforward to extend the results given for d= 1,2 to d > 2. In 
short, if p is the order of the increment used, then we are in the situation 
of Theorem A if 4 -|- 4p — 2a > d, and we are in the situation of Theorem B 
ii4: + 4p-2a<d. 

7. It may be helpful to give some intuition as to why X{t) is indepen- 
dent of the Gaussian measure B in Theorem A and the Gaussian mea- 
sure Zs in Theorem B. For simplicity consider the case d = 1 and let 
X(t) have covariance function 7(t) satisfying (1.1). Then as /i — > (cor- 
responding to a fixed-domain asymptotic regime), we have (i) Wh{t) = 
/i~"/^{X(t + h) — X{t)} converges in distribution to iV(0,2c), and (ii) 
Wh{t) is asymptotically independent of X{t). Essentially, the increments 
of B in Theorem A and Zs in Theorem B are linear combinations of terms 
of the form Wh{t), while each integrand depends only on X{t). Thus, the 
integrator and integrand are independent in the limit. In Step 6 of Sec- 
tion 5 we establish this asymptotic independence rigorously. 

8. Note that the "diagonal" {xi = X2} is explicitly excluded from the region 
of integration in Theorem B. This is in line with the definition of the 
Wiener-Ito integral given by many authors, including Dobrushin and 
Major (1979), formula (1.9), Major (1981), Theorem 8.2, Nualart (1995), 
formula (1.13) and Arcones (1994), formula (3.9), but note that in all of 
the above references the exclusion of the diagonal is not made explicit 
in the notation. See Taqqu (1979), page 77, for helpful discussion of this 
point. 

3. Numerical results. In the simulation studies described below the co- 
variance function of the underlying stationary Gaussian field was chosen to 
be of the form 7(t) = exp(— c||t||"), where a £ (0,2) and c > 0, and ||t|| is 
the usual Euclidean norm on R . We chose c = 1 when d = 1 and c = 10 
when d = 2 throughout our numerical work. The data were simulated using 
the circulant embedding approach; see, for example. Wood and Chan (1994) 
and Chan and Wood (1999) for details and further references. In all cases 
considered, it was possible to use the algorithm in its "exact" form. 

Four types of point transformations were considered: 
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1. Uniform: g{x) =^{x); 

2. Exponential: g{x) = — log{l — ^{x)}; 

3. Chi-squared with one degree of freedom: g{x) = x"^; 

4. Log-normal: (7(x;r) = exp(Tx). 

In the above $ denotes the standard normal distribution function. Two 
cases of the log-normal distribution were considered, corresponding to r = 1 
and r = 4. Note that r = 4 corresponds to an extremely nonaffine transfor- 
mation in the relevant domain and is included as an extreme case. Each of 
the above transformations preserves the fractal dimension of the underlying 
Gaussian random field [see Hall and Roy (1994)]. 

Figure 1(a) shows a realization of a Gaussian process with a = 0.1, and 
Figure l(b)-(f) shows various nonaffine transformations of this realization. 
Figure 2 shows similar nonaffine transformations of a smoother process, 
with Q = 1.0. It is clear from visual inspection of the graphs that the non- 
affine transformations do have a noticeable effect in both Figures 1 and 2, 
and in many cases the transformed processes do clearly exhibit non-Gaussian 
features. The log-normal transformation with r = 4 is particularly extreme, 
as might be expected. 

Some representative numerical results for d= 1 are displayed in Tables 1-3 
and for d = 2 in Table 4. 

Table 1 summarizes the results of a simulation study of the performance 
of three estimators, OQ^g, c^oLS ^^'^ ^GLS' °^ ^^^ fractal index a of a Gaus- 
sian process, and several transformations of this process. The notation used 
for these estimators is the same as that in Kent and Wood (1997) and Chan 
and Wood (2000): cioLg, p = or 1, is the ordinary least squares estima- 
tor of a given by (1.5), based on the log-log relationship (1.4), and using 
increment (A. 3) when p = and increment (A. 4) when p = 1; and Oqls 
is a generalized least squares estimator of a, again based on the log-log 
relationship, using increment (A. 4). For smoother processes (e.g., a = 1.9) 
the GLS estimator performs slightly better, in terms of the mean squared 
error (MSE), than the other two estimators, due to its smaller bias. For 
rougher processes (e.g., a = 1.0 and 0.1) all three estimators have very sim- 
ilar MSE with slightly higher standard deviation (SD) for the estimators 
with p= 1. One noticeable fact is that all estimators generally perform worse 
for the transformed processes. However, the deterioration is fairly mild in 
most cases, except in the case of the log-normal(4) transformation, where 
there is a substantial increase in the MSE. 

In the second study we compare the asymptotic and empirical rate of 

decrease in variance of Oqls ^^*^ '^OLS ^^ ^^® sample size n increases. The 
asymptotic rates of decrease are computed using the following results in the 
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(a) Gaussian process 




ffl|f|f|| 




W 



(b) Unjform(0,1) process 




(d) Chi-square(l) process 



uUI 


iiiMi^^ 


limm 


imMilMi 


Mm 



{e) Log-Normal(l) process 



^jJiAiUikJjii^^ 


4m 



(f) Log-Normal(4) process 



Fig. 1. Gaussian process and its point-transformed non-Gaussian processes with 
n = 1000, Q = 0.1. 
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a) Gaussian process 




(c) Exponential(l) process 




(d) Chi-square(l) process 




(e) Log-Normal(l) process 




(f) Log-Normal{4) process 




Fig. 2. Gaussian process and its point-transformed nan- Gaussian processes with 
n = 1000, a = 1.0. 



Table 1 
Comparison of the three estimators and six processes with n — 1000, m = 4, based on 100 simulations in each case 





Process 




"OLS 






•^OLS 






"gls 






a 


Bias 


SD 


MSE 


Bias 


SD 


MSE 


Bias 


SD 


MSE 




0.1 


Gaussian 


-0.021 


0.033 


0.002 


-0.021 


0.043 


0.002 


-0.021 


0.043 


0.002 






Uniform 


-0.028 


0.034 


0.002 


-0.029 


0.045 


0.003 


-0.029 


0.045 


0.003 


P 




Exp(l) 


-0.027 


0.039 


0.002 


-0.025 


0.051 


0.003 


-0.025 


0.051 


0.003 


o 

> 




2 

Xi 


-0.041 


0.041 


0.003 


-0.036 


0.054 


0.004 


-0.036 


0.054 


0.004 




Log-N(l) 


-0.032 


0.047 


0.003 


-0.029 


0.064 


0.005 


-0.029 


0.063 


0.005 


2. 




Log-N(4) 


-0.079 


0.090 


0.014 


-0.081 


0.152 


0.030 


-0.075 


0.133 


0.024 


> 
2 


1.0 


Gaussian 


-0.002 


0.041 


0.002 


0.002 


0.059 


0.003 


0.001 


0.057 


0.003 


o 




Uniform 


-0.000 


0.052 


0.003 


0.006 


0.077 


0.006 


0.005 


0.074 


0.006 


> 




Exp(l) 


-0.005 


0.055 


0.003 


-0.005 


0.076 


0.006 


-0.004 


0.074 


0.006 


i^ 




2 

Xi 


-0.008 


0.059 


0.004 


0.000 


0.079 


0.006 


0.001 


0.074 


0.005 


> 




Log-N(l) 


-0.008 


0.057 


0.003 


-0.008 


0.079 


0.006 


-0.006 


0.077 


0.006 


i 




Log-N(4) 


-0.059 


0.134 


0.021 


-0.054 


0.186 


0.038 


-0.051 


0.175 


0.033 


1.9 


Gaussian 


-0.030 


0.055 


0.004 


-0.002 


0.056 


0.003 


-0.002 


0.054 


0.003 


o 
o 




Uniform 


-0.025 


0.060 


0.004 


-0.001 


0.068 


0.005 


-0.000 


0.064 


0.004 




Exp(l) 


-0.033 


0.056 


0.004 


-0.003 


0.068 


0.005 


-0.002 


0.064 


0.004 






2 

Xi 


-0.041 


0.055 


0.005 


-0.010 


0.071 


0.005 


-0.009 


0.066 


0.004 






Log-N(l) 


-0.038 


0.054 


0.004 


-0.004 


0.072 


0.005 


-0.003 


0.067 


0.005 






Log-N(4) 


-0.068 


0.063 


0.009 


-0.013 


0.123 


0.015 


-0.012 


0.119 


0.014 





Table 2 
Comparison of the two OLS estimators and three processes with n = 1000, 2000,4000, 10,000, m — 10, based on 500 simulations in each 

case; the empirical (asymptotic) variance ratios are given when n > 2000 



Order 



Order 1 



Process 


a 


var 


2000 


4000 


10,000 


var 


2000 


4000 


10,000 


Gaussian 


0.1 


0.0^40 


0.48 


(0.50) 


0.24 


(0.25) 


0.10 


;o.io) 


0.0''60 


0.48 


;o.5o) 


0.25 


;0.25) 


0.10 


;o.io) 




0.4 


0.0^15 


0.50 


(0.50) 


0.23 


(0.25) 


0.10 


;o.io) 


0.0^21 


0.52 


;o.5o) 


0.22 


:0.25) 


0.10 


;o.io) 




0.7 


0.0^22 


0.47 


(0.50) 


0.23 


(0.25) 


0.09 


;o.io) 


0.0^32 


0.45 


;o.5o) 


0.23 


;0.25) 


0.10 


;o.io) 




1.0 


0.0^25 


0.47 


(0.50) 


0.22 


(0.25) 


0.09 


;o.io) 


0.0^37 


0.49 


;o.5o) 


0.25 


;0.25) 


0.10 


;o.io) 




1.3 


0.0^27 


0.53 


(0.50) 


0.29 


(0.25) 


0.10 


;o.io) 


0.0^36 


0.55 


;o.5o) 


0.29 


;0.25) 


0.11 


,0.10) 




1.6 


0.0^37 


0.69 


(0.57) 


0.33 


(0.33) 


0.19 


;o.i6) 


0.0^39 


0.52 


;o.5o) 


0.27 


;0.25) 


0.11 


;o.io) 




1.9 


0.0^41 


0.73 


(0.87) 


0.62 


(0.76) 


0.40 


;o.63) 


0.0''46 


0.44 


;o.5o) 


0.22 


;0.25) 


0.09 


;o.io) 


Exp(l) 


0.1 


0.0^49 


0.56 


(0.87) 


0.28 


(0.76) 


0.14 


;o.63) 


0.0^73 


0.56 


;o.87) 


0.28 


;0.76) 


0.14 


,0.63) 




0.4 


0.0^23 


0.57 


(0.57) 


0.30 


(0.33) 


0.12 


;o.i6) 


0.0^34 


0.53 


;o.57) 


0.28 


;0.33) 


0.11 


;o.i6) 




0.7 


0.0^33 


0.50 


(0.50) 


0.27 


(0.25) 


0.11 


;o.io) 


0.0^47 


0.50 


;o.5o) 


0.28 


;0.25) 


0.12 


,0.10) 




1.0 


0.0^34 


0.58 


(0.50) 


0.27 


(0.25) 


0.11 


;o.io) 


0.0^58 


0.55 


;o.5o) 


0.26 


;0.25) 


0.10 


;o.io) 




1.3 


0.0^38 


0.53 


(0.50) 


0.34 


(0.25) 


0.11 


;o.io) 


0.0^53 


0.58 


;o.5o) 


0.30 


;0.25) 


0.11 


;o.io) 




1.6 


0.0^45 


0.71 


(0.57) 


0.32 


(0.33) 


0.20 


;o.i6) 


0.0^61 


0.50 


;o.5o) 


0.28 


;0.25) 


0.10 


;o.io) 




1.9 


0.0^42 


0.74 


(0.87) 


0.63 


(0.76) 


0.41 


;o.63) 


0.0^58 


0.51 


;o.5o) 


0.25 


;0.25) 


0.11 


;o.io) 


Log-normal(4) 


0.1 


0.0^12 


0.75 


(0.87) 


0.91 


(0.76) 


0.57 


;o.63) 


0.0^33 


0.69 


;o.87) 


0.80 


;0.76) 


0.53 


;o.63) 




0.4 


o.oHi 


1.0 


(0.57) 


0.76 


(0.33) 


0.91 


;o.i6) 


0.0^21 


0.92 


;o.57) 


0.62 


;0.33) 


0.74 


;o.i6) 




0.7 


O.OHS 


0.73 


(0.50) 


0.54 


(0.25) 


0.36 


;o.io) 


0.0^36 


0.77 


;o.5o) 


0.49 


;0.25) 


0.31 


;o.io) 




1.0 


O.OH7 


0.63 


(0.50) 


0.35 


(0.25) 


0.18 


;o.io) 


0.0^37 


0.56 


;o.5o) 


0.32 


:0.25) 


0.14 


;o.io) 




1.3 


O.oHl 


0.52 


(0.50) 


0.38 


(0.25) 


0.15 


;o.io) 


0.0^29 


0.50 


;o.5o) 


0.27 


;0.25) 


0.11 


;o.io) 




1.6 


0.0^83 


0.70 


(0.57) 


0.36 


(0.33) 


0.22 


;o.i6) 


0.0^24 


0.49 


;o.5o) 


0.25 


;0.25) 


0.09 


;o.io) 




1.9 


0.0^46 


0.77 


(0.87) 


0.69 


(0.76) 


0.50 


;o.63) 


0.0^18 


0.48 


;o.5o) 


0.26 


;0.25) 


0.11 


;o.io) 



►^1 
d 

M 



O 
2 



2, 

o 

a 
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> 
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Gaussian case [see Constantine and Hall (1994) and Chan and Wood (2000)]: 

rCin-\ 0<a<3/2, 

var(aSis) ~ < C2n-^ log n, a = 3/2, 

[ C73n2"-4, 3/2 < a < 2, 



and 



var(aQL3) ~ ^4"" ' < a < 2, 



where Ci,...,C4 depend on a and m but not on n. Hence, the asymp- 
totic variance ratios for sample sizes ni < n2 and a ^ 3/2 will be equal to 
either {n2/ni)~'^ or (n2/ni)^"~^. For a non-Gaussian process in the one- 
dimensional case (d = 1), the variance formulas for both cioLS ^'^'^ '^OLS ^^^ 
the same as in the Gaussian case except when < a < 1/2. In that case the 
asymptotic variances are asymptotic to Csn"^" for both estimators. Table 2 
reports the estimated variance when n = 1000 and the estimated rate of 
decrease in variance: 

w(a„)/^^(aiooo), when n = 2000, 4000, 10,000, 

for Gaussian, exponential(l) (as an example of mildly nonafiine transfor- 
mation) and log-normal(4) (as an example of extremely nonaffine transfor- 
mation) processes. For the Gaussian case there is good agreement between 
the theoretical and the numerical results. For the non-Gaussian cases, there 
is good agreement for medium to large a (smoother process). For the ex- 
ponential(l) process, the numerical variance ratios are closer to those for 
the Gaussian process for small a (rough process). The rate of decrease in 
variance does not depend noticeably on the choice of m. 

In Table 3 we study the effect of the number of points m used in the 
regression for agLS' '^OLS ^'^'^ "^GLS' Table 3 suggests that the number of 
points used in the regression does not affect the MSE significantly for all 
three estimators. For these simulated data a choice of m = 4 would be suit- 
able. 

Next we take a closer look at the distribution of Oqls ^°^ both Gaussian 
and non-Gaussian processes. Figure 3 shows normal quantile-quantile plots 
based on 100 estimates for Gaussian, chi-square(l) and log-normal(4) pro- 
cesses. The added straight lines go through the first and third quartiles of 
these estimates and the corresponding value of the standard normal distri- 
bution. With the exception of the top right-hand corner [log-normal(4) with 
a = 0.1], all cases suggest that these empirical distributions are reasonably 
close to normal. This graphical finding is supported by the Kolmogorov- 
Smirnov goodness-of-fit test. 
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Our tentative conclusions for d= 1 are as follows: 

1. All three estimators perform fairly well under modest departures from 
normality (of the type introduced by the nonaffine transformation g). 

2. Our numerical results suggest there is no advantage in using the GLS 
estimator in the non-Gaussian case (which is not surprising, as it was 
designed for the Gaussian case). 

3. The number of points used in the regression is not critical and can be 
taken as small as 4 for the simulation data considered. 

For d = 2, the performance of the OLS estimator, based on the "square" 
increment defined in (A.8), is studied for no = (50, 50), (100, 100) and (500, 500), 
and m = 2(1)10. We denote this estimator by a„g. In particular, we com- 
pare the asymptotic and empirical rates of decrease in variance of «„„ as the 
sampling region increases. Theoretical results from Chan and Wood (2000) 
imply that for the Gaussian case, 

var(Q(„j^„^))/var(a(„2^„2)) = {n2/nif, < a < 2. 

Theorem A implies that for non-Gaussian fields, 

var(a(„,,„,))/var(«(„,,„,)) = | j^^^^^^j^.^ q < a < 1. 

Table 3 

Comparison of the three estimators for the Xi process in terms 

of MSB with n = 2000, m = 2, 4, 6, 8, 10, based on 100 

simulations in each case 





& 






m 






a 


2 


4 


6 


8 


10 


0.1 


"OLS 


0.0056 


0.0031 


0.0032 


0.0031 


0.0031 




"OLS 


0.0103 


0.0035 


0.0033 


0.0030 


0.0031 




"GLS 


0.0103 


0.0035 


0.0033 


0.0030 


0.0031 


0.3 


"OLS 


0.0054 


0.0037 


0.0041 


0.0045 


0.0048 




"OLS 


0.0111 


0.0041 


0.0040 


0.0043 


0.0045 




"GLS 


0.0111 


0.0041 


0.0038 


0.0039 


0.0040 


1.0 


"OLS 


0.0017 


0.0018 


0.0019 


0.0020 


0.0022 




"OLS 


0.0053 


0.0038 


0.0039 


0.0035 


0.0033 




"GLS 


0.0053 


0.0035 


0.0028 


0.0023 


0.0021 


1.7 


"OLS 


0.0021 


0.0025 


0.0029 


0.0032 


0.0035 




"OLS 


0.0054 


0.0033 


0.0029 


0.0032 


0.0036 




"GLS 


0.0054 


0.0030 


0.0027 


0.0024 


0.0022 


1.9 


a(0) 

"OLS 


0.0028 


0.0031 


0.0033 


0.0035 


0.0037 




"OLS 


0.0048 


0.0032 


0.0032 


0.0034 


0.0037 




"gls 


0.0048 


0.0027 


0.0024 


0.0022 


0.0023 



Table 4 
Comparison between empirical and asymptotic variance ratios among Gaussian and non-Gaussian random fields with no ~ (50,50), 

(100,100), (500,500), 771 = 4, based on 100 simulations in each case 



a 


var 


100/50 


500 


/50 


500/ 


100 


var 


100/50 


50C 


)/50 


500 


/lOO 










Gaussian field 










Uniform field 






CJ 


0.1 


0.0-= 14 


0.15 


(0.25) 


0.01 


(0.01) 


0.04 


(0.04) 


0.0^3 


0.15 


(0.87) 


0.01 


(0.63) 


0.04 


(0.72) 




> 


0.4 


0.0^15 


0.23 


(0.25) 


0.01 


(0.01) 


0.04 


(0.04) 


0.0^15 


0.21 


(0.57) 


0.01 


(0.16) 


0.05 


(0.28) 


0.7 


0.0^30 


0.24 


(0.25) 


0.01 


(0.01) 


0.03 


(0.04) 


0.0^29 


0.22 


(0.38) 


0.01 


(0.04) 


0.05 


(0.11) 


'^ 


1.0 


0.0^41 


0.20 


(0.25) 


0.01 


(0.01) 


0.04 


(0.04) 


0.0^52 


0.24 


(0.25) 


0.01 


(0.01) 


0.04 


(0.04) 


^ 
^ 


1.3 


0.0^49 


0.21 


(0.25) 


0.01 


(0.01) 


0.03 


(0.04) 


0.0^69 


0.28 


(0.25) 


0.01 


(0.01) 


0.03 


(0.04) 





1.6 


0.0^55 


0.22 


(0.25) 


0.01 


(0.01) 


0.05 


(0.04) 


O.oHl 


0.34 


(0.25) 


0.01 


(0.01) 


0.03 


(0.04) 


> 


1.9 


O.OH2 


0.26 


(0.25) 


0.01 


(0.01) 


0.02 


(0.04) 


0.0^38 


0.48 


(0.25) 


0.01 


(0.01) 


0.03 


(0.04) 


^ 










Xl field 












Log-normal(4) 


field 






> 


0.1 


0.0^15 


0.23 


(0.87) 


0.01 


(0.63) 


0.03 


(0.72) 


2.98 


0.23 


(0.87) 


0.01 


(0.63) 


0.03 


(0.72) 


1 


0.4 


0.0^20 


0.21 


(0.57) 


0.01 


(0.16) 


0.05 


(0.28) 


1.95 


0.32 


(0.57) 


0.01 


(0.16) 


0.03 


(0.28) 


0.7 


0.0^49 


0.29 


(0.38) 


0.01 


(0.04) 


0.05 


(0.11) 


3.59 


0.06 


(0.38) 


0.01 


(0.04) 


0.01 


(0.11) 






1.0 


0.0^97 


0.23 


(0.25) 


0.01 


(0.01) 


0.04 


(0.04) 


6.60 


0.23 


(0.25) 


0.01 


(0.01) 


0.04 


(0.04) 


1.3 


O.OH3 


0.25 


(0.25) 


0.01 


(0.01) 


0.04 


(0.04) 


1.36 


1.20 


(0.25) 


0.03 


(0.01) 


0.03 


(0.04) 




1.6 


O.OH6 


0.21 


(0.25) 


0.01 


(0.01) 


0.05 


(0.04) 


1.00 


0.54 


(0.25) 


0.04 


(0.01) 


0.02 


(0.04) 




1.9 


0.0^26 


0.33 


(0.25) 


0.01 


(0.01) 


0.02 


(0.04) 


2.87 


0.11 


(0.25) 


0.04 


(0.04) 


0.41 


(0.04) 
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Chi-square(l) 



Log-Normal(4) 



Quantiles of Slandard Norm 

Gaussian 



Quanliles of Standard Normal 

Chi-square(l) 



Qiianliles of Slandard Norm 

Log-Normal(4) 




^.w 



Fig. 3. Normal quantile-quantile plots for dQ^g with n = 2000, a = 0.1, 1.0, 1.9, m = 10 



Table 4 reports the estimated variance when uq = (50, 50) and the estimated 
variance ratios when ni > 712 for ni, n2 = 50, 100, 500, for the fohowing types 
of fields: Gaussian, uniform, Xi ^^id log-normal(4). For the Gaussian case 
there is good agreement between the theoretical and numerical results. For 
non-Gaussian fields there is close agreement for medium to large a. For 
small a the empirical ratios are closer to those in the Gaussian case. As in 
the one-dimensional case, the number of points used in the regression is not 
critical. In summary, the OLS estimator with p= 1 performs reasonably well 
under mild departures from Gaussianity. 

4. Some lemmas. We now present some results which are required in 
Section 5. Lemmas 4.1 and 4.2 are used to prove Lemma 4.3 which (along 
with the elementary Lemma 4.4) plays a key role in bounding remainder 
terms which arise in Steps 2-5 in Section 5. Lemma 4.5 is used in Step 6 of 
Section 5. 



Lemma 4.1. Let {Af,Bi)'^, i E Z'', where Ai = {An,. . . ,AiK)'^ , be a 
stationary Gaussian vector field and assume that {Af,Bi) has a standard 
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multivariate Gaussian distribution (i.e., with mean the zero vector and iden- 
tity covariance matrix). Define 

CTAAii - J) = , max | cov(Ajfc, Aj;)|, aBB{i - j) = I cov{Bi,Bj)\, 

(^AB [i - j) = max I cov{Aik ,Bj)\ 

k=l,...,K 

= max I cov{Bi,Ajk)\ = (JBA{i - j)- 

k=l,...,K 

Let a denote a nonnegative multi-index a= (a[l], . . . ,a[i^]) G Z , and write 
Ha{Ai) = Yli^^i Ha[k]{Aik) , where Hm, rn>0, is the Hermite polynomial of 
degree m. Then 

\E\Ha{Ai)H,„{Bi)H4Aj)H„,(Bj)]\ 



r\{aQ-ry.{aQ-r)\{m-aQ-\-r)\' 

r=max(ao— m.O) 



where qq = \a\ = J2k=i "^ W ■ 

Proof. From the diagram formula for moments (see Remark 4.1) we 
have, under the assumptions of the lemma, 

E[Ha{A,)Hra{B,)Ha{Aj)Hrn{Bj)] 

(4.1) 

ia[l]l---a[K]\mlf ^ 
= -j Z^ ^ii,ti (« - j) • • • %,t,(^ - J), 

where q = m + oq, the /'s are associated with the components of (^J", Bi)^ , 
the fs are associated with the components of {Aj,Bj)^ and the summa- 
tion is over all indices /i, ti, . . . , Iq, tq G {1, . . . , X + 1} such that there are 
precisely a[k\ ^-indices and t-indices equal to A; = 1, . . . , K, and m Z-indices 
and t-indices equal to K -\- 1. In the above, ar^s{i — j) = coY{Air,Ajs) for 
1 < r, s<K, ar,K+iii - j) = cov{Air,Bj) for l<r <K, aK+i,s{'>' ~ J) = 
cov{Bi,Ajs) and aK+i,K+iii - j) = cov{Bi,Bj). 

Consider a typical product crij^^kiii — j) " '(^lq,kq{i — j)- If this consists 
of r pairings of components of Ai with components of Aj , then there must 
be ao — r pairings of components of Ai with Bj , UQ — r pairings of components 
of Aj with Bi and m — qq -\- r pairings of Bi with Bj , where necessarily 
max(oo — m, 0) < r < ag. Therefore, 



< CTAAii - jY<TAB{i - Jr-'(TBA{i - 3r-'<TBB{i " J 
<c7AA{i-jY<rAB{i-j?'^'''-'^ 



q,r^q ^ 

ao-r^^/- A\ao-r -. (a ■\m-ao+r 
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since < dBBii — j) < 1 and (TAsii — j) = c^BAi} — j)- But the number of 
terms in the sum in (4.1) with precisely r pairings of components of Ai with 
components of Aj is bounded above by 

/ ao! W (m + ao)! 



\nf=iafc7 VJ^K^o -r)!(ao-r)!(m-ao + r)! 
Therefore, since q = m + ao, we have the bound 



(4.2) 






< 



nf=i a,\ J ,=_^_,o) ^'(«o - ^)'(«o - ^y-^-^ - «o + ^y- ' 



and the lemma follows after application of the bound (4.2) to the right-hand 
side of (4.1). D 

Remark 4.1. For details of the diagram formula for moments of prod- 
ucts of Hermite polynomials in Gaussian variables, see, for example, Taqqu (1977), 
Major (1981) and Arcones (1994). Note that the expectation on the left-hand 
side of (4.1) reduces to the expression on the right-hand side of (4.1) because 
the components of {Af,Bi)'^ are independent Gaussian variables for each i. 
Note also that the factor 2~^ which appears in Major's (1981) version of the 
formula does not appear in (4.1) because we have employed the convention 
that, for each k, Ik and tk in o"Zj.,tj. {i — j) are such that Ik is always associated 
with an z-index and tk is always associated with a j-index. 

Lemma 4.2. Let /:R^R denote a function with compact support 
whose qth derivative f^'^' is continuous on R. Write Hm{x) for the mth 
Hermite polynomial and (j){x) for the standard normal density, and let Cm = 
S^oofi^)^rn{x)(t>{x)dx denote the mth coefficient in the expansion of f in 
Hermite polynomials. Then J2m=o'^m+g/''^^- < °°- 

Proof. By assumption, f^'^' is continuous on R and has compact sup- 
port. Therefore, f^'^' has an expansion in Hermite polynomials of the form 

J2m=o Cm X Hm{x)/ml, which is L^-convergent in the sense that Z^m=o('^m )^/'^' < °^- 
But repeated integration by parts using the identity f^^(j){x)Hm{x)dx = 

—(piy) X Hm-i{y) for rTT, = 1, 2, . . . shows that Cm = Cm+q, which proves the 
lemma. D 

Lemmas 4.1 and 4.2 are used to prove the following result. 
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Lemma 4.3. Let (yl|"^^,s|"^)^, ieZ'^, n = l,2,..., be a sequence of 

stationary Gaussian vector fields, where A-^ = {Ai^ ,. . . ,Af^) is a zero- 
in) 
mean Gaussian vector whose dimension does not depend on n, and Bi ~ 

A^(0, 1). Suppose that (i) for each n Bf is independent of Af (but not 
necessarily independent of Aj when i ^ j) and (ii) the smallest eigenvalue 
of coy{A-^ ) is bounded away from as n — > co. Let 7r{Ai' ) be a polynomial 
of degree q in the components of Af^' such that E[tt{Ai'')] = for all n. 
Suppose that 

aJi(z-j)= sup |cov(4"\4f)|<C{l + |.-j|r-2p-2 

k,l=l,...,K 

and 

4t(^-J)= sup |cov(4^\i?f)|<Cn""/(2'^){l + |z-j|r-^-\ 

k=l,...,K 

where p is a nonnegative integer and C > and a G (0, 2) are constants 
independent of n. Let /i:R ^ R be expressible as a sum, of the form h = 
hi + /i2, where hi is a polynomial and /i2 is a function of compact support 
whose qth derivative is square integrable over R. Then for d=l,2, 

{0{n~^), ifp = 0and{2-a)>d, 

0{n~^logn), ifp = 0and{2-a) = d, 

0(n(°-2)/'^), tfp = 0and{2-a)<d, 

0{n-% ^fp>l. 

Proof. Assumption (ii), combined with the assumption that the ele- 
ments of cov(^^" ) are bounded above by C, imphes that we may without 

loss of generality assume that A^ is a standard multivariate normal vector 
for each n. Then, using multi-index notation again, we may write the poly- 
nomial ■k{Ai' ) as a sum of the form '^^CaLLa{Al' ), where a ranges over a 
finite set of multi-indices in Z and Ca S R- Therefore, the result will follow 
for a general polynomial vr of degree q if we can prove that it holds for each 

product HaiAl"'') = Ofc -f^a[fc] ( Afc ) o^ degree at most q. 

Let the Hermite polynomial expansion of h (which is convergent in the 
L^ sense) be given by h{x) = Y^m=o^mHm{x)/m\. Using Lemma 4.1, we 
obtain 



var| n 






OO t 

(")\ Y^ "™ tr ^r('*)^ 



var n-i 5: //.(4"^) ^ ^if„(5^ 



1 -^-^ in v -"^j 



<n-E 
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m=0 

oo 



ml 



J2 \E[Ha{Af^)Ha{Af^)Hr^{Bf^)Hrn{Bf^)] 



m=0 






(ao!)2 E ^^ 



m=0 



X E 

r=max(ao— m.O) 
ao 



r!(ao — r)\{aQ — r)\{m — ao + r)\ 



E {r!(ao -r)!(ao-r)!(m-ao + r)!} ^ 

r=max(ao— m,0) 



xn ^E^il(^-J-)VAi(^-j)^('^°-'^) 



-2 Y- ^(-) 

Using the elementary result [see Chan and Wood (2000), page 364] that 

O(n-i), iip>d 



n 



Ea + 



^-3\ 



0{n~^L{n)), ifp = d, 



where we can take L{n) = logn, and omitting some straightforward details, 
we find that the assumed bounds for a)^J^ and a^^ imply that 



n~ E -ii(^-^-)vir](.-jf(--^) 



-2 Y- ^(«), 

ro(n-i). 



ifp = 0, {2-a)>d, 
0(n~^logn), if p = 0, (2 — a) = d, 
0(n(°-2)/'^), ifp = 0, (2-a)<d, 

lO(n-i), ifp>l, 

where d = 1,2. Note that the above statement is valid for each integer m>l 
and each integer r satisfying max(0, oq — m) < r < ao- Also, using Lemma 4.2, 
it is straightforward to check that 

^ ao 



m=0 



E {r!(ao-r)!(ao -r)!(m-ao + r)!} ^ 

r=irLax(ao— "1,0) 



< OO. 



Finally, we put these results together and the proof is complete. D 

The following result is elementary but is used repeatedly, and so is stated 
explicitly for convenience. 
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Lemma 4.4. Let {X^^n'-k = l,...,n;n> 1} be an arbitrary triangular 
array of random variables such that sup^^j^ „sup„>ii?|Xfc^„| <C< oo. 
Then n~^ J2k=i ^k,n = ^^(l) as n ^ oo. 

Proof. Note that E\n~^J2k^k,n\ < C and then use the Markov in- 
equahty. 

D 

The Prohorov metric p and Ky Fan metric, here denoted k, are de- 
fined as follows [see Dudley (1989)]. Let X and Y be random elements 
of a metric space (S*, dist), with laws P and Q, respectively, defined on 
the Borel sigma field of (5, dist). Then p{P, Q) = inf{e > : P{A) < Q{A^) + 
e for all Borelsets A}, where A^ = {y £ S : dist(x, y) < £ for some x £ A} and 
k{X,Y) = inf{e > 0: P[dist(X,y) > e] < e}. Note that, by Theorem 11.3.5 
of Dudley (1989), we have 

(4.3) p{P,Q)<K{X,Y). 



Lemma 4.5. Let X and Y be random elements with laws P and Q, 
respectively, defined on the Borel sigma field of a metric space (S", dist), and 
let p denote the Prohorov metric. Then p{P,Q) < {£'[dist(X, y)^]}^'^. 



Proof. Chebyshev's inequality yields 



P 



dist(x,y) > {^[dist(x,y)^]}'/^ < {E[dist{x,Yy]} 



2ni/3 



N2ni/3 



and so the result follows from (4.3). D 

In all applications of this result given below, 5 = R^ and dist is the usual 
Euclidean metric. 

5. Proofs. We now prove the theorems stated in Section 2. The struc- 
ture of the proof of each theorem is very similar, and, in fact, the proofs 
given for Steps 1-5 cover both theorems. The only substantial difference 
between Theorems A and B is in the limit distribution which arises in Step 
6. Throughout the proof, we will use the multi-index notation specified in 
the Appendix, on the understanding that d = 1 or 2. 

We first introduce some notation that will be used throughout this section. 
Write 



,,.„...{.(^)_.(A)} 
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and 

(5.1) y^u = n-^^''^Y.-]x(^)=J:a]W.,, 

where X{t) is the underlying Gaussian field (see Section 2). Note that the 
last equality is a consequence of the fact that J2j ^^ — ^ (^se the Appendix) . 
Define 



awwii-j) = sup\cov {Wik,Wji)\, 

k,l 

<ywx{i- j) =sup|cov{VFifc,X(j7no)}|, 

k 

c^Yvii-j) =sup|cov(yi„,yj„)|, 

u,v 

cTYxii-j) = sup|cov{yi„,X(j7no)}|. 

u 

Then condition (^1)4 implies the existence of a constant C independent 
of i, j and n such that 

(Tww{i-j)<c{i+\i-j\r-^ 

(5.2) 

awx{i-j)<Cn~"^'^'''\l + \i-j\r-\ 

cTYYii-j)<C{l + \i-j\r-'P-^ 
(5.3) 

^yx(i-j)<Cn-°/(2'^)(l + |i-j|)°-P-\ 

where p is the order of the increment a on which Y is based. See Kent 
and Wood (1997) and Chan and Wood (2000), Lemma 3.1, for justification 
of (5.2) and (5.3). 

Proof of Theorem A. The proof is broken into a number of steps. 
The Ti and Tij referred to below are defined in the course of the proof; each 
of these quantities is Op(l) and in some cases of smaller order. 

Step 1. Show that it is sufficient to prove the theorem for those g which 
satisfy (.42) and (.43) and have compact support. 

Step 2. Show that 

^1(0 - an) = To + n-^/^^-i^Ti + n-^/'^Ts + Op(n-=^"/(2d) ^ ^-i)^ 
where 

(5.4) Gi = n-i^[g«{X(i/no)}]' 

iein 

and To, Ti and T2 are defined in (5.9) via (5.6)-(5.8). 
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Step 3. Show that 

To = Too + n-^/'^Toi + Op(n-3a/(2d)) 

rO(n"i/2)^ ifp = 0andd = 2, 
\ 0(n~^), if p > 1 and/or d = 1, 

where Too and Toi are defined in (5.16) and (5.14), respectively, and show 
that 



^01 ^ X! ^«Mo,m'^o,i« 



V 

•. u=l 



as n — > oo, where /zo,u and ro,iu are defined below (2.9). 

Step 4. Show that Ti = n'^'/^^'^^Tii + Op(n-"/'^), where Tn is defined 
in (5.17), and show that as n^ oo, 

m \ 



V 

V u=l 



Tn ^ ^ T«Mo,l'^o,2n 



X / b(^){X(t)}5(=^){X(t)} + [g^^\x{t)}f]dt, 

where ro,2u is defined below (2.9). 

Step 5. Show that T2 = T21 + T22 + Op(n-"/(2'^) + n-^'^), where T21 
and T22 are defined in (5.19) and (5.20), respectively, and show that, as 

n -^ 00, 

T2i^ff:T„^o>o,3«') / gW{x{t)}g^^\x{t)}dt 



and 



T22 - f E ^«/^o>o,4«') / [g^^\X{t)}fdt, 
and ro,3u and tq^^u are defined below (2.9). 

Step 6. Establish convergence in distribution of (n^'^Too, Gi), where Too 
is defined in (5.16) and Gi is defined in (5.4). 
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Proof of Step 1. Condition (.41)4 implies (1.1) when d=l and 
(2.3) when d = 2. In each case, we may use Kolmogorov's lemma [see, e.g., 
Rogers and Williams (1994), page 59] to establish that X{t) has a continuous 
version on [0, 1] . Consequently, for each e > there exists a C, depending 
on £ and the distribution of X, such that 



(5.5) P 



sup \X{t)\ >C 



<e. 



For given g, let gc denote a function with compact support such that g{t) = 
gcit) ^or all ||t|| < C, and let dc denote the estimator of a that would have 
been obtained if gc{X{t)} rather than g{X{t)} had been observed. It follows 
from (5.5) that 

P[a^ ac for some n] < e. 

As a consequence, if the theorem is true for all functions of compact support 
which satisfy assumptions (.42) and (.43), then it is also true for each g 
which satisfies (.42) and (.43), whether or not g has compact support. This 
argument can be established rigorously using probability metrics (cf. the 
argument given in Step 6). We omit the details. D 

For the remainder of the proof we shall assume that g has compact support 
[in addition to satisfying (.42) and (.43)]. 

Proof of Step 2. By Taylor's theorem 

Z] ^j'S'i+j = Z] «i (5i+i -9i) ( since Z «i = ^ ) 

= E«"fE^-™/^''^(^0-^nffr+n-^"/'^(4!)-iW^^#) 

J \r=\ / 

r=\ 

where, for r = 1, 2, 3, 



] 3 



gf is g^"^' evaluated at X(i/no), r = 1,2,3, and from Taylor's theorem, 

9ij -9i [t>j- 
ably chosen. 



g^f = gf [%X{(i + j)/no} + (1 - ej)X{i/nQ)], where each Oj G [0, 1] is suit 
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Then 

JGI„ V r=l / 

V iGJ„ / 

where 

(5.6) ri„ = 2/i-in-i ^ Mii„M2i«, 

(5.7) f2u = fi-^n~^ Y. (^4« + 2Mh„M3,„) 



JG^^n 



and 



5„ = 2(Mh,M4» + M2iuM^iu) + n-''l^^'^\Ml^ + 2M2iuMuu) 

+ 2n--/'^M3,. + M4.. + n-3-/(2'^)ML. 

Since by Step 1 we are assuming that the g^"^' are bounded for 1 < r < 4, it 
is a straightforward (if tedious) matter to check that E\Siu\ < oo. Therefore, 
since the Siu are identically distributed for each 1 <u <m, it follows from 
Lemma 4.4 that 

^-3a/(2d)^-l J2 S,u = Op(n-3a/(2d)) 

«GX,i 

for each u. Moreover, 

2 



n 



«GX„ jgX„ \ j / 



■ n 

JGXn 



— ^J'uG\ + fiuTou, 
where Gi is defined in (5.4) and 



(1)t2 



(5.8) fou = n~'Y.(^'u'Yl-l){gl'^} 

iGX„ 
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So rfl'^Zu = Gili-u + -R«/^u, where 

^« = ro„ + n-^/^^-^^fi, + n-°/'^r2„ + Op(n-3"/(2'^)). 
Since each T^^ is bounded in probabihty for k = 1,2,3, it follows that 

m 

Gi{a - an) =Gi^ L„log(Z„/^„) 

u=l 
m 

= GiY.Lu{\ogGi+\og{l + RjGi)} 

u=l 
m 

= Y,KRu{l + Op{Ru/G{)] 

u=l 
m 

w=l 

= To + n-^/^^'^^Ti + n-^/'^Ts + ©^(n-^^/^^rf) ^ „-i)^ 
where 

m 

(5.9) Tk = Y.Lufku, k = 0,1,2, 

u=l 

and we have used the fact that Ru = Op{n~^''^ + n~"' ), so that R^ = 
Op{n~^ + n~'^°' ). The order statement for R^ follows from Steps 3~6. D 

We now introduce some notation which is needed in Steps 3-5. Recall 

(n) 

the definition of Wij given at the beginning of Section 5. Writing Q = 

cov(X{i/no),W^j), define C^"^ = {Cj''\-mJ <j< mJ) and W^ = {Wij, -mJ < 

j < mJ), and let V^ denote the covariance matrix of Wj. Note that Vfy 
does not depend on i because of the stationarity of Wi, but that the dis- 
tribution of Wi does depend on n; this dependence on n has been sup- 
pressed for notational convenience. Define b^^' = {b"' ,—mJ < j < mJ) by 
5(n) ^ ^^/{2d)^y(j^)yi(-{n)_ ^^^^ ^^^^ j^-"/(2d) ^^. hf^Wij is the projection of 

X{i/n) onto the span of Wij, —mJ < j < mJ . Let b^^' = lim„^oo ^ denote 
the limit of b^^' which necessarily exists under assumption (.41)4 . 
We may write 

X{^/no) = n-"/(2'^) (llhf^W,^ + (1 - n^/'^b^-^'^V^h^-^ hof'x,. 
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where Xi = X(i/no), Xi ~ A^(0,7o) is independent of Wi, 70 = 7(0) is the 
variance of X{i/nQ) and V^^ and b^"' are as defined above. Then 

X{i/no) =X, + n~^l(^'^^u + rT'^l'^hi + n-''''l'^h^ 

where 

j 

n-^»/%, = (1 - n-"/'^6(")^vi")6(«)/7o)'/'l. - X, - n-^/%, 
and 

Proof of Step 3. Writing g\'^' for g'^^\X{i/no)) as before, and g^' 
for g^^'{Xi), we obtain 

(9i ) -[9i ) +'in 'doig^ g^ 

(5.10) 

where i?ij is a remainder term which can be determined exphcitly. 

We now study the contribution of each of the four terms on the right- 
hand side of (5.10). First, note that SQiRu can be expressed as a finite sum 
of terms, each of which can be expressed as a bounded function multiphed 
by a polynomial in Gaussian variables. Therefore, since the Sq^Ru {i €2„) 
are identically distributed, we may use Lemma 4.4 to show that 

(5.11) f: L„^-i L-' Y. (n-3-/(2'^)4i?i.)) = 0,(n-3"/(2'^)). 
u=i \ iei„ / 

Also, using similar arguments, 

2 



Es"^.)^o^a#^#^ + (#V} 



(5.12) =nA9^^gf^ + {gm 



+ {(E«"^.)'4-n«}{#^#^ + (#¥}, 



where 

(5.13) Tin = E 



E«"^d E4"V.. 
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It follows, after applying Lemmas 4.3 and 4.4 to the second term on the 
right-hand side of (5.12), that 

u=l \ ieJ„ \ j / / 

where 



m 



^(l)x(3),/x(2)^2 



(5.14) Toi = E Lut^Z'riu U"' E 9^9^ + (^l 

\«=1 / V i€ln / 

A similar argument, using Lemmas 4.3 and 4.4 again, shows that 

(5.15) n-"/(2'^)n-i E (E^.^^v) ^o.^f ^fff ^ = ©.(n'^), 

except when p = and d = 2, in which case we can conclude that the left- 
hand side of (5.15) is of size Op(n~^'^). Writing 

(5.16) Too = E Lun-' Y^i^r^Vl - l){g^^f 

U=l i&Tn 

and putting (5.10)-(5.16) together, we see that the first part of Step 3 is 
proved. To establish the limiting distribution of Tqi, we may use an identical 
argument to that given in Step 6 to prove that Jn{h) -^ Jo{h) [which is the 

more straightforward part of showing that Kn{h) -^KQ{h)]. To avoid dupli- 
cation, we omit the details. 

D 

Proof of Step 4. By a Taylor expansion. 

Using similar arguments to those in Step 3, in particular Lemma 4.3, we 
find that for any e > 0, 

Ti = E ^«/^« ^ {--' E f E -1^^^ f E «^^i) ^^ 

«=1 \ JGXn \ j / \ k I I 

where 

(5.17) Tn= f E^«^«'-2„)n-i E{#^#^ + (#V} 

\«=1 / JGI„ 
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and 

(5.18) T2U = e[ fe«i^«.) fe«^^i) (E^S"^^^') I' 

so the first part of Step 4 is proved. Justification of the claim concerning the 
limit distribution of Tn follows along similar lines to that for Tqi in Step 3. 

D 

Proof of Step 5. In this case 

Jl)^(3) _-(!)- (3) ,^-a/(2d)r p ^. ('^(2) ^ 2 _ .^(2) x 2 -a/(2d) . p 

9i 9i —9i3i +n 'OoiKii and [gi ) —[gi ) +n 'OoiH^i. 

Using Lemmas 4.3 and 4.4 again and writing 

(5.19) T21 = [^L./z-Va^'jn-i ^ t^gf\ 

(5.20) r22= (E^-^nV)n-i 5:(#V, 

(5.21) r3„ = ii?fea|Ty,,'| fe«i^)' ^^u = \E{y^a^Wf\ , 

we find that n~°^''^T2 = n~°^''^i^2i +^22) + Op{n~^) as required. Justification 
of the claim concerning the limit distributions of T21 and T22 follows along 
similar lines to that for Tqi in Step 3 and Tn in Step 4. D 

Proof of Step 6. Here we shah show that 

(5.22) (ni/2Too,Gi)^('a / [g^^^{X(t)}f dB{t),Gi 

where Tqq is given in (5.16), Gi is defined in (5.4), Gi is defined in (2.5) with 
r = 1, {B{t)} is the random Gaussian measure given in (2.11), which is inde- 
pendent of {X(i)}, and a is given in (2.10). Recall that no = ('^oil]) . . . j'^oi^]) 
and n = ni=i^o[^] and that we assume that condition (^4) in Section 2.1 
holds when d = 2. 

Let 7i denote the class of smooth functions with compact support. For h G 
7i, define 

In{h) = n^i/2 J2 ^ih{X{i/no)}, Jn{h) = n~^ ^ /i{X(i/no)} 

and K^{K) = (I„(/i), J„(h)), where ^i = E™=i Luif^-^Yii - 1). Note that by 
construction K„,{(g(i))2} = {n^/'^Too,Gi). Also define 

Io{h)= f h{X{t)}dB{t), Mh)= f h{X{t)}dt 

J[0,1]'* "'[0,1]'* 
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and Ko{h) = (/o(/i), Jo(/i)). Note that by construction KoKfi' )^} is equal 
to the right-hand side of (5.22). 

We wiU show that, for each h€7i, Kn{h) converges to KQ{h) in distribu- 
tion. Then, in view of (5.22) and Step 1, Step 6 will follow. 

For each positive integer r, let vr^ denote a partition of [0, 1]'' given by 
TTj. = {Aj,j e Jr}, where j;. = {j G Z'^ : < j[l] <r,l<l <d}, where Aj C 
[0, 1]"^ is defined in the following way. Write tj = r~^j. When d = 1, Aj is 
the interval of width r~^ centered at tj = r~^j + g; when d = 2, Aj is a 
square with sides of length r~^ which are parallel to the coordinate axes, 
and each Aj is centered at tj = r~^j + 2^, where 1 = (1, 1) is a 2-vector 
of I's. Given iTr, we define two functions, t*{t) and i*{i), as follows: for 
t G [0, 1] , t*(t) = tj when t ^ Aj, and for each multi- index < i < no, we 
define i*{i) = ini^) — ^o^ji where j is such that i/riQ G Aj. Define 

and K*Xh) = {I*{h),J*{h)), and write 

I^{h)=I^{h;7rr)=a f h{X{t*)}dB{t), 

Jo* = Jo(/i;vrr)= / h{X{t*)}dt 

J [0,1]'* 

and KQ{h) = (/q(/i), Jq (/i)), where {B{t)} is a random Gaussian measure of 
the form (2.11) which is independent of {X{t)}. 

Let Pn,hi Po,h, Pnh ^^^ ^^ h denote the distributions of Kn{h), KqQi), 

K*{h) and KQ{h), respectively. We need to show that Pn,h^ Po,h- We shall 
do this by showing that, given e > 0, there exists a partition Tr^ of [0,1]^^, 
depending on e and h, such that 

(5.23) p{Po,h,Po,h)<e/3 
and, for any such vr,., when n is sufficiently large, 

(5.24) p{Pl^,Pl„)<e/3; 

and for a sufficiently fine partition vr,. and n sufficiently large, 

(5.25) p{Pn,h,P:,h)<^/^, 

where p denotes the Prohorov metric. Then, when n is sufficiently large, 

p{Pn,h, Po,h) < p{Pn,h, Pn,h) + PiKh, P^h) + P(^0,h, ^0,0 < £, 
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and, since e > may be chosen arbitrarily small, Step 6 will have been 
proved. 

Proof of (5.23). Using Lemma 4.5, 

(5.26) p{Po,h,P^,h) < [E{h{h)-r,{h)Y+E{Uh) - r^{h)ff\ 

Using Fubini's theorem, the Cauchy-Schwarz inequality and the fact that {X[t)} 
and {-B(t)} are independent, we obtain 

E{h{h) - rQ{h)Y = eL f [h{x{t)}-h{x{t*)}]dB{t) 



(5.27) 

'[0,1] 
and 



a^ I E[h{X{t)} - h{X{t*)}fdt 
^fo.il<* 



(5.28) E{jQ{h) - JQ{h)Y < f E[h{X{t)} - h{X{t*)}fdt. 

J [0,1]'* 

Given h £7i and e > 0, it is clear (since {X{t)} is almost surely continuous 
and therefore uniformly continuous on [0, l]'^) that we can find a (sufficiently 
fine) partition tt^ of [0,1]'^ such that 



/ E[h{X{t)}-h{X{f)}?<U^/3y 
"'foal'* 



^3/2 
'[0, 
Then, using (5.27) and (5.28), 

E{Io{h)-I^{h)}'<'^{e/3f/^ and E{Jo{h) - J^{h)}^ < '^{e/3f/^ , 
in which case (5.23) follows from (5.26). D 

Proof of (5.24). For any Lebesgue-measurable set AC [0,1]"^, define 
(5.29) En{A) = n~'/^ Y. ^- 

Then, for a given partition vr,., 

/:(/i) =n-V2 ^ i,h{X{i*/n,)] = Y: E4A,)h{Xit,)} 

i&In j&Jr 

and 

^: W = n'^ Y HX{i*/no)} = n-i ^ c,h{X{t,)}, 

where Cj = ^{i:i/nQ G Aj} ~ n/r . In all cases covered by Theorem A a 
central limit theorem applies to H.„([0, 1]''); that is, H„([0, 1]'') ^A^(0, cr^). 
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where o^ is defined in (2.10). See Kent and Wood (1995, 1997) and Chan and 
Wood (2000) for further details. It follows from a straightforward extension 

of those proofs that, for any fixed vr^, {H„(Aj), j S Jr] -^{aB{Aj),j S J'r}, 
where B is the random Gaussian measure defined in (2.11). Moreover, it is 

an easy consequence of the definition of X that {X{tj),j G J'r}—>-{X{tj),j G 
J'r}, where the random variables {X(tj),j E J'r} are independent of the ran- 
dom variables {B{Aj), j £ Jr}. Consequently, an application of the continu- 

ous mapping theorem implies that K*[h) — > K^Qi), from which (5.24) follows 
for sufficiently large n, since the Prohorov metric p metrizes convergence in 
distribution. D 

Proof of (5.25). We will only sketch the proof of this result. It is 
sufficient to show that 



'^dt 



(5.30) lim E{lJh) - I* (h)}^ = a^ f E[h{X(t)} - h{X(t*)}] 

n^oo "'[0,1]'' 

and 

(5.31) \im EUJh) - J*(h)}^ < f E[h{X(t)} - h{X(t*)}f dt, 

because, if we choose a partition vr such that, for all n sufficiently large, 
/ E[h{X{t)} - h{X{t*)}f dt < i(e/6)3/2, 

J[0,l]d 

then (5.25) will follow from (5.30), (5.31) and Lemma 4.5. 

The proof of (5.31) is omitted as it is quite straightforward. However, 
we outline the proof of (5.30), which requires more work. For simplicity, 
we focus on the case in which h = J2m'^mHm is a polynomial; it is an easy 
matter to extend the proof to all h£7i. Now 

E{Uh) - /:(/i)}2 = n-i Yl ^fe^i(^i - h^'){hJ - h,,)] 



M 

— — -n > 6{i,j,m,m), 

, _, mim'l . '^—i 

m,m'=l i,jGl„ 



where hi = h{X{i/nQ)}, hi* = h{X{i* /uq)} and 

6{i,j,m,m') = E[i,ij{Hm{X,) - Hm{X,,)}{H^,{X,) - H^,{Xj,)}]. 
Then (5.30) is a consequence of the following: for each m ^ m' , 
(5.32) n~^ J2 5ii,J,m,m')=o{l); 
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and 

lim n~ > 6(i,j,m,m) 



(5.33) 



n — >oo ^ — ^ 

X [HmiXj) - Hm{Xj*)] 

: a^ f [H^{X{t)) - Hm{X{t*))f dt. 



'[0,1]' 

The results (5.32) and (5.33) are derived using the diagram formula (see 
the references given in Remark 4.1) combined with (5.2) and (5.3). The 
arguments are broadly similar to those used in the proof of Lemma 4.3. 
That concludes our sketched proof of (5.25). D 

Thus, Step 6 is now complete. D 

To conclude the proof of Theorem A, we use (5.22) and the continuous 
mapping theorem to show that 

GrV/^Too^f / [,«W)}]^d5(t)^.^Z, 

where Z ~ A^(0, 1) is independent of Gi and G2, and G2 is defined in (2.5) 
with r = 2. Finally, putting Steps 1-6 together, we see that the proof of 
Theorem A is complete. D 

Proof of Theorem B. The proof of Theorem B is essentially the same 
as the proof of Theorem A, except that Step 6 is different. The decomposition 
given by (5.23)-(5.25) can still be used, but the principal difference is that 
S„ in (5.29) does not obey a central limit theorem, and therefore In{h) does 
not converge to a stochastic integral with respect to the random Gaussian 
measure B{t). What actually happens is that, when d=\ and 3/2 < a < 2 
or d = 2 and 1 < a < 2, 



var n 






e =0(,l(2"-4)/d)^ 



and n'^~"" "^n"^ X^jgx ^« obeys a noncentral limit theorem of the type given 
by Theorem 6 in Arcones (1994). Then it can be shown, via the decomposi- 
tion (5.23)-(5.25), that In{h) converges to the Wiener-Ito integral specified 
in the statement of Theorem B. The particular form of the spectral mea- 
sure 5 follows in part as a consequence of Theorem 1' of Dobrushin and 
Major (1979) and in part as a consequence of the degeneracy result given in 
part (c) of Theorem 1 in Kent and Wood (1997). D 
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APPENDIX: NOTATION FOR INCREMENTS 

Further details concerning increments may be found in Kent and Wood (1997) 
and Chan and Wood (2000); we only give brief details here. 

The univariate case (d = 1) . An increment of order p > is a finite array 
of real numbers a = {aj : —J <:j<:J} such that, for all integers 0<r<p, 

(A.l) ^ faj = 0, 

j--J<j<J 

but 

(A.2) ^ f+\ / 0. 

j--J<j<J 

An example of an increment of order p = is given by 

(A. 3) a = {ao,ai}, where ao = —1 and ai = 1; 

an example of an increment of order p = 1 is given by 

(A. 4) a = {a_i,ao,ai}, where a_i = l,ao = ~2 and oi = 1. 

Note that, here and in the main text, we adopt the convention that any 
components aj which are not given explicitly are [so, in (A. 3), for example, 
we have omitted a„i = 0]. 

For an integer u> 1, the dilation a" = {a'' : —Ju < j < Ju} of an incre- 
ment a has elements defined by 

(A.5) <=l?'' 'ir^''""' 

■^ [U, otherwise, 

where —Ju < j < Ju. It follows immediately from this definition that 

ro, ifO<r<p, 

^/aj= n-^^/a,, ifr>p+l, 
J [ j 

where here and in the main text ^ • indicates summation over —Ju < j < Ju. 
Let {nj :j€'Z} be a sequence of numbers. Then if the increment a is given 
by (A. 3), it follows that 

/ , O-j Hn+j ^ Vn+u ~ Vn j 
j 

while if a is given by (A. 4), then 

^ Ojyn+j = Vn+u + Vn-u " '^Vn- 
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The multivariate case {d> 1). Let j = (j[l], . . . ,j[d\) G Z denote a multi- 
index. We say that j is nonnegative and write j > if j[l] > for each I = 
1, . . . ,d; and if k is another multi-index, we say that j < k if j[l] < k[l] for 
each /, and write j <k ii each inequality is strict. For multi- indices j and r, 
we define 

d 
(A.6) \r\=y2r[l] 



=1 



and 



(A.7) f = EM'^'^^ 

1=1 

where 0*^ = 1. 

In the d-dimensional case, an increment of order p > is a finite array 
a = {ttj : —J < j < J} satisfying (A.l) and (A. 2), but with j, J and r now 
multi- indices with j^ defined by (A.7) and each aj real as before. 

The so-called "square" increment in the case d = 2, which is an example 
of an increment of order p= 1, is given by 

(A.8) a = {a(o,o) = a(i,i) = 1, «(i,o) = a(o,i) = -1}; 

see Chan and Wood (2000). 

The dilation a" = {a^ : —Ju < j < Ju} is defined by (A. 5), but with j and 
j' now multi-indices. It follows immediately from this definition that 

where \r\ is given by (A.6). 

Note that if {yj : j G Z^} is an array of real numbers and the square in- 
crement (A.8) is used, then for any n G Z^, 

/ . O'^Vn+j = yn+{0,0) + yn+{u,u) ~ Un+iufi) ~ ?/n+{0,u)- 
j 
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